February 1, 2008 
BA-TH/99-359 



Spontaneous emission and lifetime modification 
caused by an intense electromagnetic field 

P. FACCHI and S. PASCAZIO 

Dipartimento di Fisica, Universita di Bari 
and Istituto Nazionale di Fisica Nucleare, Sezione di Bari 

1-70126 Bari, Italy 

PACS: 42.50.Hz; 42.50.Vk; 03.65.Bz 
Abstract 

We study the temporal evolution of a three- level system (such as an atom or a molecule), 
initially prepared in an excited state, bathed in a laser field tuned at the transition frequency 
of the other level. The features of the spontaneous emission are investigated and the lifetime 
of the initial state is evaluated: a Fermi "golden rule" still applies, but the on-shell matrix 
elements depend on the intensity of the laser field. In general, the lifetime is a decreasing 
function of the laser intensity. The phenomenon we discuss can be viewed as an "inverse" 
quantum Zeno effect and can be analyzed in terms of dressed states. 



1 Introduction 



The temporal behavior of quantum mechanical systems can be strongly influenced 
by the action of an external agent. A good example is the quantum Zeno effect 
[0, [|, where the quantum mechanical evolution of a given (not necessarily unstable) 
state is slowed down (or even halted) by performing a series of measurements that 
ascertain whether the system is still in its initial state. This peculiar effect is his- 
torically associated and usually ascribed to what we could call a "pulsed" quantum 
mechanical observation on the system. However, it can also be obtained by perform- 
ing a "continuous" observation of the quantum state, e.g. by means of an intense field 

II- 

Most experiment that have been performed or proposed in order to modify the 
quantum mechanical evolution law make use of oscillating systems f| |6|, [7|, || || . On 
the other hand, it would be interesting to understand whether and to which extent 
the evolution law of a bona fide "unstable" system can be changed. In order to discuss 
the evolution of genuine unstable systems one usually makes use of the Weisskopf- 



Wigner approximation [10|, that ascribes the main properties of the decay law to a 



pole located near the real axis of the complex energy plane. This yields the Fermi 



"golden rule" ||TTJ. In this paper we shall investigate the possibility that the lifetime 
of an unstable quantum system can be modified by the presence of a very intense 
electromagnetic field. We shall look at the temporal behavior of a three-level system 
(such as an atom or a molecule), where level #1 is the ground state and levels #2, 
#3 are two excited states. (See Figure ffl.) The system is initially prepared in level 
^2 and if it follows its natural evolution, it will decay to level #1. The decay will 
be (approximately) exponential and characterized by a certain lifetime, that can be 
calculated from the Fermi golden rule. But if one shines on the system an intense 
laser field, tuned at the transition frequency 3-1, the evolution can be different. This 
problem was investigated in Ref. ||, where it was found that the lifetime of the 
initial state depends on the intensity of the laser field. In the limit of an extremely 
intense field, the initial state undergoes a "continuous observation" and the decay 
should be considerably slowed down (quantum Zeno effect). The aim of this paper 
is to study this effect in more detail and discuss a new phenomenon jnj: we shall 
see that for physically sensible values of the intensity of the laser, the decay can be 
enhanced, rather than hindered. This can be viewed as an "inverse" quantum Zeno 
effect. An important role in this context will be played by the specific properties of 
the interaction Hamiltonian, in particular by the "form factor" of the interaction. 

Other authors have studied physical effects that are related to those we shall 
discuss. The features of the matrix elements of the interaction Hamiltonian were 



investigated in the context of the quantum Zeno effect by Kofman and Kurizki []T3 
who also emphasized that different quantum Zeno regimes are present. Plenio, Knight 
and Thompson discussed the quantum Zeno effect due to "continuous" measurements 
and considered several physical systems whose evolution is modified by an external 
field []14||. There is also work by Kraus on a similar subject flT5[. Finally, Zhu, Narducci 



and Scully |16| investigated the electromagnetic-induced transparency in a context 



1 



similar to that considered in this paper. In some sense, our present investigation 
"blends" these studies, by taking into account the important role played by the matrix 
elements of the interaction. This will enable us to discuss some new features of the 
evolution that have not been considered before. We shall look at this phenomenon 
from several perspectives, by first solving the time-dependent Schrodinger equation, 
then looking at the spectrum of the emitted photons and finally constructing the 
dressed (Fano) states. 

Our analysis will be performed within the Weisskopf-Wigner approximation and 
no deviations at short |], [17], |18| and long |1| times will be considered. The features 
of the quantum mechanical evolution are summarized in and have already been 
discussed within a quantum field theoretical framework [[H], 22 , 23, , where several 
subtle effects have to be properly taken into account. 

This paper is organized as follows: in Section 2 we introduce the 3-level system 
bathed in the laser field. Its temporal evolution is studied in Section 3. The spectrum 
of the photons emitted during the evolution is evaluated in Section 4. Section 5 
contains a discussion in terms of dressed states, Section 6 an analysis of the influence 
of additional levels on the lifetime and Section 7 some concluding remarks. 



2 Preliminaries and definitions 

We consider the Hamiltonian (h = c = 1)@: 
H = H + H int 

= co 1 2} (2 1 + Q 1 3} (3 1 + J2 WfcoLafcA + E (<^a4a l 1 )( 2 l + (j)* hX a kX \2) (1 1 

fc,A fc,A 

+ E (**a°LJ1><3| + $* feA a feA |3)(l|) , (2.1) 

fc,A 

where the first two terms are the free Hamiltonian of the 3-level atom (whose states 
\i) (i = 1,2,3) have energies E x = 0, u = E 2 - E 1 > 0, £l = E 3 - E 1 > 0), the 
third term is the free Hamiltonian of the EM field and the last two terms describe the 
1^2 and 1 «-> 3 transitions in the rotating wave approximation, respectively. (See 
Figure [!].) States |2) and |3) are chosen so that no transition between them is possible 
(e.g., because of selection rules). The matrix elements of the interaction Hamiltonian 
read 



^fcA 



PC Ps 



J d x e e "*> ' Mx) •» = TE^S ' d x e «*» ' ^<* : 



(2.2) 

where — e is the electron charge, eo the vacuum permittivity, V the volume of the 
box, uj = \k\, ek\ the photon polarization and j & the transition current of the radi- 
ating system. For example, in the case of an electron in an external field, we have 
3 a = ipfOcipi where ip\ and ip{ are the wavefunctions of the initial and final state, 
respectively, and ct is the vector of Dirac matrices. For the sake of generality we 
are using relativistic matrix elements, but our analysis can also be performed with 
nonrelativistic ones j & = iplpipi/m ei where p/m e is the electron velocity. 
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Figure 1: Level configuration 



We shall concentrate our attention on a 3-level system bathed in a continuous 
laser beam, whose photons have momentum ko (|&o| — ^o) an d polarization Ao, and 
assume, throughout this paper, that 



^fcoAo 



0, (2.3) 



i.e., the laser does not interact with state |2). Let the laser be in a coherent state 
|ao) with a very large average number No = \ao\ 2 of fco-photons in volume V [we will 
eventually consider the thermodynamical limit; see Eq. ( p.l 7| )1. In the picture defined 
by the unitary operator 

T(t) = exp (a*e inot a ko x ~ a e- iUot a{ oXo ) , (2.4) 
the Hamiltonian (|2.1| ) reads 

H(t) = THT^ + ifTl = H+ ($ fcoAo a*e^|l>(3| + <^ oAo a e-^|3)(I|) • ( 2 - 5 ) 

In this picture, the ko mode is initially in the vacuum state [ES] and by noting that 
for N > 1 

0fcx|#(*)|3; fcA )| = V^I^qAoI > \(l;l k x\H(t)\3;0 kX )\ = |$ fcA |, (2.6) 
the Hamiltonian (|2.5| ) becomes 

H ~ ^ o |2)(2| + fi o |3)(3| + ^^aLa fc A + E'(^A4A|l)(2|+0L«fcA|2)(l|) 

fc,A fc,A 

+ ($ feoAo a*e^*|l)(3| + $* oAo a e-^|3)(l|) , (2.7) 

where a prime means that the summation does not include (ko, Ao) [due to hypothesis 
Q2.3Q 1. In the above equations and henceforth, the vector \i;rik\) represents a state 
in which the atom is in state \i) and the electromagnetic field in a state with n k \ 
(fe, A)-photons. We shall analyze the behavior of the system under the action of a 
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continuous laser beam of high intensity. Under these conditions, level configurations 
similar to that of Figure p] give rise to the phenomenon of induced transparency |]26 



for laser beams of sufficiently high intensities. Our interest, however, will be focused 
on unstable initial states: we shall study the temporal behavior of level #2 when the 
system is shined by a continuous laser of intensity comparable to those used to obtain 
induced transparency. 

Notice that in Eq. ( |2.7|) the spontaneous decay 3 — >• 1 has been neglected with 

respect to the stimulated transition, because of the large factor \fJSfo ^> 1 in Eq. 
(|2.6|) . However, since our interest is primarily in the first step of this process, namely 
the decay 2 — > 1, these smaller, later effects (of the order of 1/N ) do not change our 
conclusions. 
The operator 

M=\2)(2\+J2'4 x a kX , (2.8) 

fe,A 

satisfies 

[H,Jtf\ = 0, (2.9) 

which implies the conservation of the total number of photons plus the atomic exci- 
tation (Tamm-Dancoff approximation P7|). The Hilbert space splits therefore into 
sectors that are invariant under the action of the Hamiltonian: in our case, the sys- 
tem evolves in the subspace labelled by the eigenvalue AT = 1 and the analysis can 
be restricted to this sector pf[ . 

3 Temporal evolution 

We will study the temporal evolution by solving the time-dependent Schrodinger 
equation 

ij t \m) = H(t)\m), (3.i) 

where the states of the total system in the sector M = 1 read 

= x(t)\2; 0) + J2'ykx(t)\l; l fcA ) + J2'zkx(t)e~ inot \3; l kX ) (3.2) 

fc.A fc,A 

and are normalized: 

(mm)) = w*)i 2 + EW*)i a + eWwi 2 = 1. w (3.3) 

fc,A fc,A 



By inserting ( |3.2j ) in ( |3.1| ) one obtains the equations of motion 
i x(t) = uox(t) + E <Pl\yk\{t), 

k,X 

iyk\{t) = <Pkxx(t) + uJ k y kX {t) + a* ® koX[) z kX (t), 

iz kX (t) = a $ koXo y kX (t) + uj k z kX (t), (3.4) 
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where a dot denotes time derivative. At time t = we prepare our system in the 
state 

|V(0)) = |2; 0) ^ x(0) = 1, y kX (0) = 0, ^ fcA (0) = 0, (3.5) 
which is an eigenstate of the free Hamiltonian 

ifc#(0)> = #o|2;0> =^ |2;0). (3.6) 

Incidentally, we stress that the choice of the initial state is different from that of Ref. 
||, where the 3-level atom is initially in the ground state (#1) and a Rabi oscillation 
to level #2, provoked by an rf-field, is inhibited by a pulsed laser, resonating between 
levels #1 and #3, that performs the "observation" of level #1. In our case, the 
atom is initially in level #2, so that it can spontaneously decay to level #1, and it is 
"continuously observed" by a continuous laser at the 1-3 frequency 0: As soon as 
the system has decayed to level #1, the (intense) laser provokes the 1 — > 3 transition. 
(The irreversibility inherent in the act of observation is eventually brought in by the 
spontaneous decay of level #3.) This brings us conceptually closer to the seminal 
formulation [T], of quantum Zeno effect. 

By Laplace transforming the system of differential equations ( |3.4|) and incorpo- 
rating the initial condition ( ^.5| ) we get the algebraic system 

isx(s) = uj x(s) + J2 ( l ) *k\yk\(s) +i, 

k,\ 

is Vk\{s) = (f)kxx(s) + u k y kX (s) + a*^ ko \ Zk\{s) , 

isz kX {s) = a ®* koXo y k \(s) +uj k z kX (s), (3.7) 

where 

f(s) = / dt e~ st f(t), fit) = — da e ts f(s), (3.8) 

JO 2 717 Jb 

the Bromwich path B being a vertical line Res =constant in the half plane of conver- 
gence of the Laplace transform. (Very similar equations of motion can be obtained 
by assuming that the external (laser) field is initially in a number state N , with N 
very large [pL2 ] . See also the discussion in Section |5|.) It is straightforward to obtain 



%( s ) = : 1 \ ? (3-9) 



ykx(s) = - , ... , x 2 - m x{s), (3.10) 



-i4>k\(s + iv k ) 
{s + tu k ) 2 + B 2 

^fc) 2 + 



with 



and where 



^ (s + %<jj k y + ±? 2 



fe,A 



B 2 = N \$ koXo \ 2 (3.13) 
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is proportional to the intensity of the laser field and can be viewed as the "strength" 
of the observation performed by the laser beam on level #2 ||. See the paragraph 
following Eq. (|3.6|) . Note that the coupling B is related to the Rabi frequency by the 
simple relation B = fi Rabi /2. 

In the continuum limit (V — > oo), the matrix elements scale as follows 



v Vu 2 
hm 



V-^oo (27r) 3 



J2 J ^|0 fc A| 2 = ^oxV), (3-14) 



where Q is the solid angle. The (dimensionless) function x( a; ) an d coupling constant 
g have the following general properties, discussed in Appendix A: 



cu 2 ^ 1 if uj < A 
g 2 = a(u Q /A) 2j+lTl , (3.16) 



ifw>>A , (3.15) 



where j is the total angular momentum of the photon emitted in the 2 — > l transition, 
=F represent electric and magnetic transitions, respectively, /3(> l) is a constant, a 
the fine structure constant and A a natural cutoff (of the order of the inverse size of 
the emitting system, e.g. the Bohr radius for an atom), which determines the range 



of the atomic or molecular form factor 29 



In order to scale the quantity B, we take the limit of very large cavity, by keeping 
the density of f2 -photons in the cavity constant: 

V — > oo, N — > oo, with — = n = const (3-17) 



and obtain from ( |3.l3D 



B 2 = n V\$ koXo \ 2 = (2nyn \<p Xo (k )\ 2 , (3.18) 

where (p = $V l / 2 /(27r) 3//2 is the scaled matrix element of the 1-3 transition. As we 
shall see, in order to affect significantly the lifetime of level #2, we shall need a 
high value of B, namely, a laser beam of high intensity. It is therefore interesting to 
consider a 1-3 transition of the dipole type, in which case the above formula reads 

B 2 = 2rran \e* koXo -x l3 \ 2 n 0j (3.19) 

where X13 is the dipole matrix element. 



3.1 Laser off 

Let us first look at the case B = 0. The laser is off and we expect to recover the 
well-known physics of the spontaneous emission a two-level system prepared in an 
excited state and coupled to the vacuum of the radiation field. In this case, Q(0, s) 
is nothing but the self-energy function 

Q(s) = Q(0, s) = |0fcA| 2 ^— , (3.20) 

fc,A S + %UJk 
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Figure 2: Cut and pole in the s-plane (B = 0) and convergence circle for the ex- 
pansion of Q(s) around s = —iujQ. I and II are the first and second Riemann sheets, 
respectively. The pole is on the second Riemann sheet, at a distance 0(g 2 ) from 
-iu . 



which becomes, in the continuum limit, 



roc "^(lo\ 

Q(s) = g 2 uj q(s) = -ig 2 u / duj -. 

J o u — IS 



(3.21) 



where x is defined in (|3.14f) . The function x(s) in Eq. ( p.9[) (with B = 0) has a 
logarithmic branch cut, extending from to —ioo, and no singularities on the first 
Riemann sheet (physical sheet) [p2] . On the other hand, it has a simple pole on the 
second Riemann sheet, that is the solution of the equation 



s + iu + g 2 u Q qii(s) = 0, 



(3.22) 



where 

q u (s) = q(se- 2 ™) = q(s) + 2n X 2 (ts) (3.23) 

is the determination of q(s) on the second Riemann sheet. We note that g 2 q(s) is 
0(g 2 ), so that the pole can be found perturbatively: by expanding qn(s) around — iu 
we get a power series, whose radius of convergence is R c = u because of the branch 
point at the origin. The circle of convergence lies half on the first Riemann sheet and 
half on the second sheet (Figure |^). The pole is well inside the convergence circle, 
because |s po i e + iu>o\ ~ g 2 u <ti R c , and we can write 

Spoie = -iuj -g 2 uj q n (-iuj -0 + ) + O(g' i ) = -iu -g 2 u q(-iuj + + ) + O(g 4 ), (3.24) 

because qn(s) is the analytic continuation of q(s) below the branch cut. By using the 
formula 

lim = P-T mS(x\ (3.25) 

e^0+ X ± is X 
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one gets from ( |3.21| ) 

1 



q(—ir] + + ) = —i dux 2 {u) 



u; — 77 — i0 + 

= irx 2 {r))6{r))-iPj Q du X 2 (u;)-j-^ (3.26) 

and by setting 

Spoie = -iuo + iAE - ^, (3.27) 

one obtains 



7 = 2irg 2 u oX 2 M + 0{g% AE = g 2 u P H du -^M. + 0(<? 4 ), (3.28 

JO id — LUq 



which are the Fermi "golden rule" and the second order correction to the energy of 
level #2. 

The Weisskopf-Wigner approximation [|nj consists in neglecting all branch cut 
contributions and approximating the self-energy function with a constant (its value 
in the pole), that is 

1 1 1 , 

X(S) = ; — - ~ ; — = , (3.29) 

s + iuq + Q(s) s + iu Q + Qn[s po]ie ) s - s po i e 



where in the last equality we used the pole equation ( |3.22| ). This yields a purely expo- 
nential behavior, x(t) = exp(s po i e t), without short-time (and long-time) corrections. 
As is well known, the latter are all contained in the neglected branch cut contribution. 



3.2 Laser on 

We turn now our attention to the situation with the laser switched on (B 7^ 0) and 
tuned at the 1-3 transition frequency Q . The self energy function Q(B, s) in ( |3.12| ) 
depends on B and can be written in terms of the self energy function Q(s) in absence 
of laser field [Eq. ( |3.2(J| )], by making use of the following remarkable property: 

Q(B, s) = W |0 fc A| 2 ( —■ 1 ■ - R + * -J = \ m + iB) + Q(s - iB)] . 
2^ V s + iuOk + iB s + iuJk — iB J 2 

(3.30) 

Notice, incidentally, that in the continuum limit (V — > 00), due to the above formula, 
Q(B,s) scales just like Q(s). The position of the pole s po i e (and as a consequence 
the lifetime te = 7 _1 = — l/2Res po i e ) depends on the value of B. There are now two 
branch cuts in the complex s plane, due to the two terms in ( |3.30| ). They lie over the 
imaginary axis, along (— zoo, — iB] and (— ioo, +iB]. 
The pole satisfies the equation 

s + ilu + Q(B,s) = 0, (3.31) 



S 




Figure 3: Cuts and pole in the s-plane (B ^ 0) and convergence circle for the ex- 
pansion of Q(B, s) around s = —iujq. I , II and III are the first, second and third 
Riemann sheets, respectively, (a) B < ojq. (b) B > ojq. In both cases, the pole is at 
a distance 0(g 2 ) from — ilu . 

where Q(B,s) is of order g 2 , as before, and can again be expanded in power series 
around s = —iouo, in order to find the pole perturbatively. However, this time one has 
to choose the right determination of the function Q(B, s). Two cases are mathemati- 
cally possible: a) The branch point —%B is situated above — iu , so that — iu lies on 
both cuts. See Figure |3|(a); b) The branch point —iB is situated below — iu , so that 
— iu lies only on the upper branch cut. See Figure |](b). We notice that, although 
mathematically conceivable, the latter case (B > ujq) cannot be tackled within our 
approximations, for a number of additional effects would then have to be considered: 
multi-photon processes would take place, the other atomic levels would start to play 
an important role and our approach (3-level atom in the rotating wave approxima- 
tion) would no longer be valid. We therefore restrict our attention to values of B 
that are high (of the same order of magnitude as those utilized in electromagnetic 
induced transparency), but not extremely high, so that our starting approximations 
still apply. 

In case a), i.e. for B < Uq, the pole is on the third Riemann sheet (under both 
cuts) and the power series converges in a circle lying half on the first and half on the 
third Riemann sheet, within a convergence radius R c = uj — B , which decreases as 
B increases [Figure 0(a)]. For the sake of completeness we also notice that in case 
b), i.e. for B > co> , the pole would be on the second Riemann sheet (under the upper 
cut only) and the power series would converge in a circle lying half on the first and 
half on the second Riemann sheet, within a convergence radius R c = B — ujq, which 
increases with B [Figure 0(b)]. 
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In either cases we can write, for |s po i c + iwo\ < R c = \B — u \, 

-iw - ~ {QH(lo + B)+ 0+] + Q[-i(u - B) + + ]} + 0(g*) 
-iu - ]-g 2 uo Q {q[-i{uoo + B) + 0+] + g[-i(w - 5) + 0+]} + 



'polo 



Equation ( p.32| ) enables us to analyze the temporal behavior of state #2. 



(3.32) 



3.3 Decay rate vs B 

We write, as in (|3.27|) , 



Spoie = -iuo + iAE(B) 



l{B) 



(3.33) 



Substituting ( |3.26| ) into ( |3.32| ) and taking the real part, one obtains the following 
expression for the decay rate 



j(B)=irg 2 u x 2 {^ + B)+x 2 {^-B)d{uu Q -B) +0(g 4 ). 



On the other hand, by (|3.28[ ), one can write 

X> + B) + X > - B)8(u -B) 



l{B) =7 



+ 0(g 4 



(3.34) 



(3.35) 



This is the central result of this paper and involves no approximations: Equation 
( |3.35|) expresses the "new" lifetime r y(B)~ 1 , when the system is bathed in an intense 
laser field B, in terms of the "ordinary" lifetime 7 -1 , when there is no laser field. 
By taking into account the general behavior ( 3.15|) of the matrix elements x 2 ( w ) an d 
substituting into ( |3.35| ), one gets to 0(g A ) 



1(B) 



1 
2 



1 + 



B\ 2 ^ 



B_\ 2 ^ 



8(uj - B) 



[B < A) (3.36) 



where =F refers to 1-2 transitions of electric and magnetic type, respectively. Observe 
that, since A ~ inverse Bohr radius, only the case B < uj A is the physically 
relevant one \T2\. The decay rate is profoundly modified by the presence of the laser 
field. Its behavior is shown in Figure f| for a few values of j. In general, for j > 1 (1-2 
transitions of electric quadrupole, magnetic dipole or higher), the decay rate j(B) 
increases with B, so that the lifetime 7(5) _1 decreases as B is increased. If one looks 
at B as the strength of the "observation" performed by the laser beam on level ^2 
[[J, one can view this phenomenon as an "inverse" quantum Zeno effect, for decay is 
enhanced (rather than suppressed) by observation. 

As we shall see in Sections |] and |5], the emitted photons have different frequencies 
[for they correspond to decay onto different dressed (Fano) states]. By selecting the 
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i{B)h 



i{B)h 




photon wavelength (i.e., by means of filters), one could therefore also measure the 
different contributions to the inverse lifetime in (|3.35|) - (|3.36|) . We shall come back to 
this point later. 



As already emphasized, Eq. ( 3.36 ) is valid for B <C A. In the opposite (unphysical) 
case B > A, by Q3TT5D and (|0|), one gets to 0(g 4 ) 

1(B) * \ ^| oc (B/A)-P. (B » A) (3.37) 

This result is similar to that obtained in Ref. [[|. If such high values of B were 
experimentally obtainable, the decay would be considerably hindered and B could 
be properly viewed as the "strength" of the observation performed by the laser field 
on level #2 (quantum Zeno effect). However, in such a case, many additional effects 
would have to be considered and our analysis should be modified in order to take 
them into account. A similar remark was made by Kofman and Kurizki in a different 



context 13 



A final remark is now in order. If one would use the approximation ( |3.29|) in Eq. 
( |3.30| ), in order to evaluate the new lifetime, i.e. if one sets Q(s) = Q(s po i e ) =const, 
one would obtain Q(B,s) = Q(s) = Q(s po \ e ), i.e. no .B-dependence. Therefore, the 
effect we are discussing is ultimately due to the nonexponential contributions arising 
from the cut. In particular, viewed from the perspective of the time domain, this 
effect is ascribable to the quadratic short-time behavior of the 2 — > 1 decay. 



3.4 Estimates 

We saw in the previous subsection that the ratio B/u Q is the relevant quantity in the 
evaluation of the modified lifetime. Let us therefore try to get a rough feeling for the 
magnitude of the relevant physical parameters. In order to affect significantly the 
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lifetime of level #2, we have to look at rather large values of B: for instance at 1-3 
transition of the electric dipole type. In such a case, Eq. (|3.19|) applies: 



B 2 = 27raQ \e* koXo ■ x 13 \ 2 n . (3.38) 
Considering the angle average 

(l4oA -*i 3 r> = ^i3| 2 (3.39) 
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and remembering that the decay rate is 

4 



r 13 = ^a\x l3 \ 2 Q 3 , (3.40) 



we obtain 



5 2 = |n ^, (3.41) 



which, reinserting c's and fr's, reads 

b 2 = |n M]o^r 13 = (n nn )j^(hr 13 ), (3.42) 

where Xl = 2ttc/VLq. The quantity B 2 has dimensions of squared energy and is given 
by the product of the energy of the laser field contained in the volume A|/167r 2 times 
the energy spread of the 1 — 3 transition Qq. Therefore B depends on both laser and 
atomic system. Observe that uqX 3 l is the number of laser photons contained in the 
volume A|. 

In terms of laser power P and laser spot area A, Eq. (|3.42|) reads 



P A 3 PA 3 

b 2 = ^x^t 2 ^ = 132 ^r^) eV2 > ( 3 - 43 ) 

where P is expressed in Watt, Xl in fim, A in /im 2 and for in eV. In Eq. (|3.43|) 
the quantity B is expressed in suitable units and can be easily compared to u [the 
ratio B/ujo being the relevant quantity in Eq. ( |3.36| )1. For laser intensities that are 



routinely used in the study of electromagnetic induced transparency, the effect should 
be experimentally observable. For a quick comparison remember that B is just half 
the Rabi frequency of the resonant transition 1 — 3 [see paragraph following Eq. 
TT3D1. 



4 Photon spectrum 

It is interesting to look at the spectrum of the emitted photons. It is easy 
to check that, in the Weisskopf-Wigner approximation, the survival probability 
|(^(0)| , ?/'(t))| 2 = |x(t)| 2 decreases exponentially with time. The standard way to 
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dP/dcu 



dP B I duj 





Figure 5: The spectrum (|4.5 ) of the emitted photons. The height of the Lorentzians 
is proportional to the matrix element x 2 ( w ) (dashed line). We chose an electric 
quadrupole transition, with j = 2 and 7 = lO -1 ^, and used arbitrary units on the 
vertical axis, a) B=0; b) B = lu /5; note that from 7(5) = (28/25)7. 



obtain this result is to neglect the cut contribution in the complex s plane, or equiv- 
alently, to substitute in (|3.9|) the pole determination of the self-energy function: 



1 1 1 / ^ 

nj^l g \ — ~ — M ]7 

s + iu + Q(B, s) s + iu} + Q(B,s pole ) s-s po i e ' 



from which one gets 

( - l( B ) \ 
x(t) = exp(sp i e t) = exp I -iu t —tj , (4.2) 

where Uq = ujq — AE(B). In this approximation, for any value of B, the spectrum 
of the emitted photons is Lorentzian. The proof is straightforward and is given in 
Appendix B. One finds that, for B = 0, the probability to emit a photon in the range 
(lj, oj + du) reads 

dP B =o = g 2 uJ X 2 (u)f L (uj - lu ; i)du } (4.3) 

where 

h(<»\i)= * m - (4-4) 

OJ + 7/4 

On the other hand, when B 7^ one gets: 

dP B = 9 2 ujoX 2 (^ [f L (uj -u -B; j(B)) + f L (u - oo + B; j(B))] duo. (4.5) 

The emission probability is given by the sum of two Lorentzians, centered in ujq ± B. 
We see that the emission probability of a photon of frequency ujq + B (lDq — B) 
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increases (decreases) with B (Figure |^). The linewidths are modified according to 
Eq. Q3.36Q . When B reaches the "threshold" value ujq, only the photon of higher 
frequency (u>o + B) is emitted (with increasing probability vs B). 

Photons of different frequencies are therefore emitted with different rates. We 
shall understand better the features of the emission in the next section, by looking 
at the dressed states of the system. 



5 Dressed states and links with induced trans- 
parency 

It is useful and interesting to look at our results from a different perspective, by 
analyzing the modifications of the energy levels provoked by the laser field. For 
simplicity, since the average number N of fc -photons in the total volume V can be 
considered very large, we shall perform our analysis in terms of number (rather than 
coherent) states of the electromagnetic field. In this limit, 

(l;0 fcA ,iV |# int |3;0 fcA ,iV -i> = //Vo$ feoAo > (1; Ux, N -l\H- mt \3; fcA , N -l) = $ fcA , 

(5.1) 

V(fe, A) 7^ (k , X ). [This is equivalent to (|2.6|) .] In the above equation and henceforth, 
the vector |i;nfc A ,M ) represents an atom in state with n^x (k, A)-photons and 
Mq laser photons. 

In the above approximation, the Hamiltonian ( |2.1| ) becomes 



H ~ ^ |2)(2| + fi |3)(3|+^^aLa fc A + E , (^A4A|l)(2| + ^ fc A|2)(l| 

k,X k,X 

+ (^ fco Ao4 A |l)(3| + $* feo Ao a fcoA |3}<l|) , (5.2) 

where a prime means that the summation does not include (k , Xq) [due to hypothesis 
31)]. Besides (|2.8|), there is now another conserved quantity: indeed the operator 



A/"o=|3)(3| + 4 oAo a fcoAo (5.3) 

satisfies 

[H,M ] = [Mo,M] = 0- (5.4) 

In this case, the system evolves in the subspace labelled by the two eigenvalues Af = 1 
and Ao = iV , whose states read 

|V>(t)) =x(t)|2;0,iV )+E , 2 / fc A(t)|l;l fc A,iV ) + ^ fcA (t)|3;l fcA ,iVo-l). (5.5) 

fc,A fc,A 



By using the Hamiltonian (|5.2j ) and the states ( f>.5\ ) and identifying A^o with A^o = 
|ao| 2 j the Schrodinger equation yields again the equations of motion (|3.4|) , obtained 
by assuming a coherent state for the laser mode. Our analysis is therefore indepen- 
dent of the statistics of the driving field, provided it is sufficiently intense, and the 
(convenient) use of number states is completely justified. 
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Energy conservation implies that if there are two emitted photons with different 
energies (as we saw in the previous section), there are two levels of different energies to 
which the atom can decay. This can be seen by considering the laser-dressed (Fano) 
atomic states |3D |. The shift of the dressed states can be obtained directly from 
the structure of the Hamiltonian ( |5.2|) . In the sector Ao = iV , the operator Af is 
proportional to the unit operator, the constant of proportionality being its eigenvalue. 
Hence one can write the Hamiltonian in the following form 

H = H — Q Af + Q N , (5.6) 

which, by the setting Ei + N Q = 0, reads 

H = H + H int 

= W 1 2) (2 1 + ^kCtkxCtkX + £' (^a4 a |1><2| + ^ A a feA |2><l|) 
fc,A fc,A 

+ ($ fc oAo4oA |l)(3| + ^ Ao«feoAo|3)(l|) • (5.7) 

On the other hand, in the sector H_\r^f with Af = 1 and J\f = N , the last term 
becomes 

(^ Ao4oAo|l)(3| + ^ Ao^oAo|3)(l|) = (<& feo Aa V^ol 1) (3 1 + $l X V^) (1 1) ■ 

(5.8) 

Let us diagonalize this operator, i.e. let us look for two non-interacting states |+) and 
|— ) which are linear combinations of the old states |1) and |3). To this end we write 



ID = ^ (1+) + *»!-> 
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72 



+) 



-)) 



with |+) e |— ) orthonormal: 

(+!+> = H-> 



Plugging ( |579D into 
H m t = 



i, (+R = o. 

, the interaction term becomes 



fc,A 

*( 

+B 



V2 



b k\\ 



->(2| + %a fcA |2><- 



>k\ AS t 



V2 



e a 



-)(2| + 



'fcA 



J(3 



V2 

(l+>+e"|->) ((+ 



OfeA|2)(- 



-iSi 



(-1) +h.c. 



where we have set $k \oV~No = Be %13 . Rearranging the last term 



B 



(l+X+l -!-)(-! + e l5 \-)(+\-e- l5 \+)(-\) 



(5.9) 



(5.10) 



(5.11) 
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+)(+!- |->H+ e" i5 |+)( 



AS I 



-)( + 



Bcos(P - a) (l+X+l - l-X-l) + iBsin(/3 - a) (e i<5 |-)(+ 



-><- 



(5.12) 



and setting a — (3 the two states |+) and |— ) decouple and one gets 



h 



mt 



E' 

fc,A 

+ u 

+B\+)( 



k\ f 
a kX 



+)(2| + ^a fcA |2)(+ 



V2 



kX i5 t 
e a k\ 



)(2| 
^|-}(- 



v/2 
ft* 

7T 



5 a fc A|2)(- 



Therefore we can write 



Ho + -f/int 

where the primed free and interaction Hamiltonians read respectively 
H' = LU \2)(2\ + B\+)(+\-B\-)(-\+J2'u k ai x a kX , 



(5.13) 
(5.14) 



fc.A 



hL 



E' 

fc,A 



^fcA f 



y/2 



4aI+)(2| 



yfcA 
y/2 



OfcA|2)(+| 



^ai x |->(2| + ^a h x|2>< 



(5.15) 



and we set 5 = 0. We see that the laser dresses the states |1) and |3), which (if one 
includes the Qo photon) are degenerate [with energy E — 0, due to the choice of the 
zero of energy: see line after (5J3)]. The dressed states |+) and |— ) have energies +B 
and —B and interact with state |2) with a coupling 4>k\/V2- Since 2B = f^Rabi these 
are the well-known Autler-Townes doublet []3j]] . 

Therefore, by applying the Fermi golden rule, the decay rates into the dressed 
states read 

, 2, rr, , 2, , r?> 

(5.16) 



7+ 



2 x 2 K-g) 

27T5- 



2 X 2 K + ^) 

7_ = 2^ 



and the total decay rate of state |2) is given by their sum 

7 = 7+ + 7-, 



(5.17) 



which yields ( |3.34|) . One sees why there is a threshold at B = ojq\ For B < u , 
the energies of both dressed states |±) is lower than that of the initial state \2) 
[Figure 0(a)] . The decay rate 7_ increases with B, whereas 7+ decreases with B; 
their sum 7 increases with B. These two decays (and their lifetimes) could be easily 
distinguished by selecting the frequencies of the emitted photons, e.g. by means of 
filters. 

We also notice, for completeness, that when B > u , the energy of the dressed state 
|+) is larger than that of state |2) and this decay channel disappears [Figure ||(b)]. 
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\2>- 

l+>- 
ll>- 

h>- 




(a) 



- «o 

,+B 
. 
,-B 




(b) 



Figure 6: Shift of the dressed states |+) and |— ) vs B. (a) For B < loq there are two 
decay channels, with 7_ > 7+. (b) For B > ujq level |+) is above level |2) and only 
the 7_ decay channel remains. 



As repeatedly emphasized, this situation is unphysical and would require a different 
analysis, for additional effects would play an important role. 

Finally, let us emphasize that if state |2) were below state |1), our system would 
become a three-level system in a ladder configuration, and the shift of the dressed 
states would give rise to electromagnetically induced transparency p6| . The situation 



we consider and the laser power required to bring these effects to light are therefore 
similar to those used in induced transparency. 

6 Influence of other levels 

Let us now see how our results are modified by the presence of off-resonant levels. 
To this end we generalize the three-level Hamiltonian ( f2.7| ) by including other off- 
resonant levels I j) (j = 4, • • • , N) in our analysis: 

TV 

H = ^ |2)(2| + ^ n,|j><j| + ^^at A a fc A + (^> fcA aL|l><2| + ^* A a fcA |2> <1|) 

j=3 k,\ k,\ 

N 

+ E (%a* e^\l)(j\ + $*« e-^|j)(l|) , (6.1) 
i=3 

where $j = $j.fc A are the matrix elements of the 1 <-> j transitions and Qj = Ej — E% 
the energy of level \j) [in particular, $3,fc a — 

$ fcoAo and Q 3 = n in Eq. (gTfy]. 
By a calculation similar to that used in Section |3] one gets again the expression 
(|3.9D for the Laplace trasform of the survival amplitude, with the new self-energy 
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function modified by the presence of other levels 



>fcA 



N 

J=3 s+i8j+iu! k 



(6.2) 



M s + za;/- + X! 

where = |$j| 2 /|$ 3 | 2 and <Sj = fi, - fi 3 . 

The denominator of Q(-B, s) is now a polynomial of order N — 1 (when N = 3 one 
reobtains Eq. (|3.12|) with a quadratic polynomial). Hence the new Q(B,s) in (|6.2|) 
has N — 1 branching points and the property ( |3.30|) is generalized to 



Q(B, s) = c+Q(s + ia + ) + c^Q(s + io_ 



N 
J=4 



1(7 



J I ' 



(6.3) 



where {— ia + , — ia-, —i<Jj} (j = 4, ■ • • , N) are the branching points, i.e. the zeroes of 
the denominator of Q(B,s). In this case one has to solve an algebraic equation of 
(N — l)th order, whose zeroes do not have in general an analytical expression. We 
seek a perturbative solution in B. It is lengthy, but straightforward, to obtain up to 
second order in B 



a± = ±B 



r>2 v^TV fj 
2^=4 25, 



c ± = l/2 T BT,Uit 



B 2 E N h 



}2fj 



(6.4) 



From the above equations we see that the presence of off-resonant levels modifies the 
energies a± = ±B of the two dressed states by a shift of order B 2 and creates N — 3 
new dressed states with energies 5j ~ Qj — fi 3 , whose contribution to the self-energy 
function is of order B 2 . 

By starting with the self-energy function (|6.3|) and looking for the location of the 
pole one obtains instead of Eq. (|3.36|) the following expression for the modified decay 
rate 



Tmany (-^) 



7 



w . 

N 

e=4 



OA 



c 1 +c+ 1-^ 0(w o -o- + ) 



^0 

6»(^ - cr<?) 



(6.5) 



where k — 2j =p 1. 

By substituting the expressions ( |6.4|) for the zeroes and the coefficients, valid up 
to second order in B, into Eq. (|6.5Q one gets 



Tmany (-^J 



7 S 1 + «~2 + 772 J2 ft 



uJ 



o e=4 



K 



UJ 




( w 


uj 2 


k 


5j) 


U 





d2 N 

7 (S) + 7 — Y,h 



o e=4 



K 



UJ 


w 2 \ 


( w 


UJ 2 


h 


5j) 







9{uq - 5 1 
6(u - S t 



+ 0(5 3 ), 
(6.6) 
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where ^y(B) is the decay rate ( |3.36| ) evaluated in the three-level approximation. 

The above general expression can be evaluated in practical cases of interest. For 
instance, by assuming that the off-resonant levels are well separated from the three 
main levels, that is by assuming 5i = Qi — Q% > u , all dressed states other than |±) 

because their energies are larger than the energy uj of level 



do not enter in Eq. 
1 2), and Eq. (|6.6p reads 



7many(5) ~ 7 



(6.7) 



where 



B* = B 



N 

Yh— 

k Mi. 



(6.1 



This is the correction sought: the effect of sufficiently off-resonant levels, 5g > luq, 
modifies the decay rate (|3.36 ), calculated in the three-level approximation, simply by 
changing B into B* . Observe that fe is a rapidly decreasing function of I (polynomial 
fall-off in atomic systems). Notice also that B* > B, hence the presence of the other 
levels enhances the effect discussed in Section ^. 



7 Concluding remarks 

We have studied the evolution of an unstable system under the action of an external 
(laser) field. The dynamical evolution of level #2 (initial state) is modified by the 
laser field, tuned at the transition frequency 1-3. For physically sensible values of 
the parameters, the decay of level #2 is faster when the laser is present. Equations 
( |3 . 35| )- ( pT3^ ) (valid to 4th order in the coupling constant) express the new lifetime as a 
function of the "natural" one and other parameters characterizing the physical system. 
The initial state decays to the laser-dressed states with different lifetimes. We have 
obtained Eq. ( |3.35| ) in 3 different ways, deriving the Fermi golden rule from the time- 
dependent Schrodinger equation, by making use of Laplace transforms, as in Section 
|3~3| , or starting from the dressed states, as in ( |5. 16|) - (pTT7| ) , or as a consequence of a 



normalization condition, as in ( B.13 ). We also computed, in Section 6, the corrections 
due to off- resonant levels. We emphasize that, since we always work in the Weisskopf- 
Wigner approximation, the conceptual problems related to state preparation ]22j and 
deviation from exponential behavior p0| , |Tj, [ll| were not considered. 

In which sense is the phenomenon discussed in this paper an "inverse" quantum 
Zeno effect? If the situation B ^> A were experimentally attainable, then decay would 
be hindered and one could reasonably speak of a quantum Zeno effect provoked by 
the "continuous" observation performed on the system by the laser beam. On the 
other hand, when 5 C A, one can still think in terms of a "continuous gaze" of the 
laser on the system, but this enhances (rather than hinder) decay. One should also 
notice that the inclusion of the spontaneous decay of level #3 in the Hamiltonian 
would not change our conclusion (up to order Ti 3 /B). The interpretation in 
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terms of an "inverse" quantum Zeno effect is appealing and enables one to look at 
the problem from a different perspective. 

Acknowledgments We thank E. Mihokova, G. Scamarcio and L.S. Schulman for 
useful discussions. 



Appendix A 

We discuss here some general properties of the matrix elements and derive Eqs. Q3.15Q - 
fl3.16|) of the text. An exhaustive analysis of some general features of the matrix 



elements can be found in [p9 |, but we will focus here on the behavior at small and 



large values of to. The matrix elements (|2.2j) of the 1-2 interaction Hamiltonian read 

d 3 xe~ lkx e* kX .j l2 (x), (A.i; 



12-Ka 



Vu 



where a = e 2 /4neo is the fine-structure constant. If the wavelength of the radiation is 
large compared to the size a of the system (i.e. u <C A = a -1 ) the main contribution 
to the integral ( |A.1|) comes from small values of r = \x\ (lot < 1). Expanding the 
exponential (k = nu) 

( iuS] 2 

exp(—ik ■ x) = exp(— icon ■ x) = 1 — iuj{n ■ x) H — — (n • x) 2 + ... (A. 2) 

and integrating term by term one obtains the asymptotic series 



oo 



where 

q { :l= [ -^fjd 3 xe* nX .j 12 (x)(n-xy (A.4) 

(e^x = e n \ depends only on the direction of k). Notice that we explicitly wrote every 
cu-dependence and that does not depend on oj. Observe that g^ ^ corresponds to 
electric-dipole transitions El, g^ 1 ) to electric quadrupole E2 and/or magnetic dipole 
transitions Ml, and so on. Hence s = j — A, where A = (A = 1) stands for magnetic 
(electric) transition Mj (Ej). Since the dominant contribution to the integral in ( |A.4|) 
comes from a region of size a and the current ji 2 is essentially uq/o 2 , we get 

£l<xuJoa s+1 , lu = E 2 -E 1 . (A.5) 



If ooa 1 the dominant term in the series (|A.3| ) is the first nonvanishing one, namely 



(Pkx ~ <?2^ r oc (uj Q a)(uja) r , (A. 6) 
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for some s = r. In the continuum limit one gets 



^fcA 



DC 



doo 9 ojqX (w), (A.7) 



where we have defined 



sWV) = Hm £ j rf^ 1 fcA | 2 , (Ai 



as in (|3.14|) . From (|A.6|) we obtain 



r=0 

and therefore 



,2, . „ ,2/ 



WW ~ |£/ dOlffSl 2 ) ^ 2r+1 « K^) 2r+ Vo (^) 2r+1 • (A.10) 

Remembering that 2r + 1 = 2j — 2 A + 1 = 2j =p 1, we obtain the first equation in 
( P15Q and Eq. ( g3§ . 

On the other hand, if the wavelength is much smaller than a (i.e. ^> A), we 
first rewrite (|A.l ) in the following form 



kX= ^j d3X e ' WnXe ^ • = V / ^'W^ll), (All) 



where 



JnA,i2(^[|) = / d 2 x ± e* nX ■ j 12 (x) (A.12) 

and x = xyn + ajj,. According to the Riemann-Lesbegue lemma, the integral in ( |A.11| ) 
vanishes in the uj —>■ oo limit. In particular, if Jt*a,i2 (^h ) is N times different iable, 
integrating by parts we get 



Vuj (ioj) 



N 



1 2ixa If i,,„d , . , k s 

J dx e-^ x —j nX , 12 {x) (A.13) 



and we can write 

4> kX = o{LO- N - l l 2 ), («>A) (A.14) 

which yields the large uj behavior of the second equation in ( |3.15|) . It goes without 
saying that if Jtt.a,i2(^||) is an analytic function, then ip^x — > more rapidly than any 
power. The second equation in ([3.15 ) is therefore a conservative estimate. 
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Appendix B 



In this Appendix we shall analyze the spectrum of the emitted photons. We start by 
substituting (|4.1|) into fl3.10|) and (|3.11|) , to obtain 



Vkx{s) 



Zk\(S 



-i<j>k\(s + iuJ k ) 1 

(s + iu k ) 2 + B 2 S- Spole 
\/^0 $ fcoAo^A 1 



(s + iuj k ) 2 + B 2 s - Spoic' 
Closing the Bromwich path with a semicircle in the half plane Res < 0, we get 



(B.l) 
(B.2) 



Vk\(t) 



— f ds e ts y kX {s), 

— f ds e ts z k \(s), 
m Jr 



(B-3) 



which can be evaluated by summing over the integrand residues. The quantity 
\ykx(t)\ 2 (\zkx{t)\ 2 ) represents the probability that, at time t, the transition 2 — > 1 
(2 — > 1 — > 3) has taken place. When t — > oo, the contribution of s po i e (that has a 
finite negative real part) is exponentially damped. This leaves only the contributions 
of the poles in —i{uJf. ± B). 

We look first at the case B = (laser off). One gets (z k x = 0, Vt) 



\ykx(+oo) 



?feA 



(uj k - io ) 2 + 7 2 /4 



(B.4) 



and, in the continuum limit ( 3.14|) , the probability to emit a photon in the frequency 
range (u, u + duj) reads 



dP B =o = g 2 u x 2 {^)h{^ - ^o; j)du), 
where Jl is the Lorentzian profile 

h(u;j) = J— 1 2/ ■ 
ur + 7^/4 



(B.5) 



(B.6) 



This is Eq. ( ^4.3| ) of the text. The quantity P must be normalized to unity: imposing 
this condition one gets the Fermi golden rule ( p.28|) . 

On the other hand, when B ^ 0, the total emission probability is given by the 
sum 

|y fcA (oo)| 2 + |^ fcA (oo)| 2 (B.7) 

and it is straightforward to derive the following expressions (uk = uj^ — uoq and we 
write for simplicity "f(B) = 7) 



\ykx{oo) 



>kX\ 



B< 



7 



75 



V k + J r) cos 2 (Bt) + B sin 2 (Bt) + sin(2fii) 
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>k\\ 



B< 



1 



v 2 k + -r) sin 2 (Bt) + B 2 cos 2 (Bt) - — sin(2£t) 



which yield 



(B.8) 



\ykx{00)\ + \z kx {oo) 



; fcA| 



fa + *i 



B 1 



7 
4 



(B.9) 



Therefore, in the continuum limit, we can write 



dP B = g 2 uJ X 2 {u) 



[uo-uoof + ^ + B^ 



[lu - lu - B) 2 + ^ (oj-u + B) 2 + 



duj. 



(B.10) 



This formula can be rewritten in the following form 

dP B = g 2 u x 2 (u})^ [f L (u - u) Q - B; 7) + f L (u -u> + B; 7)] du. 



(B.11) 



This is Eq. ( |4.5| ) of the text. We see that the emission probability is the sum of two 
Lorentzians, centered in lj — B and loq + B and weighted by g 2 u)oX 2 ( UJ )- This result 
is in agreement with that obtained in Refs. |16| , |32]| . Incidentally, we notice that the 
value Q3.34Q of j(B) can be readily estimated by imposing the normalization of the 
emission probability 

/r°° I 
dP B = J o g 2 oJoX 2 ^)^[fL^-oJ -B- n ) + f L (u;-6j + B^)}du = l. (B.12) 

Performing the integration one obtains (7 ujq, hence one can integrate over the 
whole real axis and take x 2 ( w ) equal to its value on each Lorentzian peak) 



1 r i 2lT 

I = / dP B « -g 2 uj [x 2 (iO + B) + X 2 (i0 - B)\ 



(B.13) 



which yields Eq. ( |3.34| ) of the text. 
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